External-Compression Supersonic Inlet Design Code 

Dr. John W. Slater, NASA Glenn Research Center / Inlet and Nozzle Branch (RTE) 



A computer code named SUPIN has been developed to perform aerodynamic design and analysis 
of external-compression, supersonic inlets. The baseline set of inlets include axisymmetric pitot, 
two-dimensional single-duct, axisymmetric outward-turning, and two-dimensional bifurcated-duct 
inlets. The aerodynamic methods are based on low-fidelity analytical and numerical procedures. 
The geometric methods are based on planar geometry elements. SUPIN has three modes of 
operation: 1) generate the inlet geometry from a explicit set of geometry information, 2) size and 
design the inlet geometry and analyze the aerodynamic performance, and 3) compute the 
aerodynamic performance of a specified inlet geometry. The aerodynamic performance quantities 
includes inlet flow rates, total pressure recovery, and drag. The geometry output from SUPIN 
includes inlet dimensions, cross-sectional areas, coordinates of planar profiles, and surface grids 
suitable for input to grid generators for analysis by computational fluid dynamics (CFD) methods. 
The input data file for SUPIN and the output file from SUPIN are text (ASCII) files. The surface 
grid files are output as formatted Plot3D or stereolithography (STL) files. SUPIN executes in batch 
mode and is available as a Microsoft Windows executable and Fortran 95 source code with a 
makefile for Linux. 
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Background 


Goal : Develop computational tools to perform aerodynamic design and analysis of 
supersonic inlets to determine inlet geometry and performance. 

Some key points: 

o Supersonic inlet design methods have been developed over the last 60+ years. 

o Various computational tools have been developed over the decades (IPAC, InletMOC, 
Lerclnlet, LAPIN) based on analytic, empirical, and computational methods. 

o New computational frameworks (Java, Matlab, etc...) and new inlet concepts (stream- 
traced, advanced flow control, etc...) have prompted us to re-visit our supersonic inlet 
design tools and to develop new tools. 

o The tools should perform low-fidelity analysis while also providing geometry for higher- 
fidelity, computational fluid dynamics (CFD) analysis. 

o This has lead to the Inlet Tools effort of the Supersonic Project. 

o Small-scale effort: the work discussed in this presentation has involved one researcher 
(John Slater) at V 2 full-time employee (FTE) for about 3 years. 
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Four Aspects to the Inlet Tools Effort 



1. Inlet Tools Electronic Database 

• Internet-based eRoom (access is controlled by member account 
with userid and password). 

• Database stores files submitted by members. 

• Database contains PDF files of public-domain papers, reports, 
and other literature. 

• Location to store data, source code, input and output files. 

2. Inlet Design and Analysis Documentation 

• Documents fundamental aerodynamic design and analysis methods suitable for supersonic inlets. 

• Provide background for methods coded into the computational tools. 

• Document the inlet geometry model and computational tools. 

• The document is in the form of a MS Word document. 

3. Inlet Geometry Model 

• Describes supersonic inlet geometry in terms of basic parts and key geometric factors. 

• Recognizes that some inlet geometry can be constructed from simple, planar entities. 

• Provides framework for constructing surfaces for visualization and CFD analysis. 

4. Computational Tools 

• Implements the inlet geometry model and design and analysis methods into a computer program. 

• Various frameworks can be used (Fortran 95, Java, Python, C++, Matlab, etc...) 


Aspects of the Inlet Tools Effort 

1 . Electronic Database 

2. Design and Analysis Document 

3. Geometry Model 

4. Computational Tools 
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SUPIN (SUPersonic INIet) Design Code 


• Incorporates the supersonic inlet geometry model and design and analysis methods. 

• Currently capable of design and analysis of external-compression, supersonic inlets (Machl .6-2.0). 

• Written in Fortran 95 on Windows laptop. 

• Text-based (ASCII) input data file with blocks and lines of input. 

• Text-based (ASCII) output file containing inlet coordinates and performance properties. 

• Creates surface and planar CFD grid files (ASCII Plot3d and STL formats). 

• Single code with an Windows executable or source code and makefile for Unix / Linux systems. 

• Modes of inlet design and analysis: 

° K mode = 1 ■ Generate the inlet geometry from an explicit set of geometry factors and inputs. 

° K m0 de= 2 - Perform design operations to size and generate the geometry of the inlet based on 
specified factors and inputs. 

° K m0 de= 3- Perform aerodynamic analysis of an inlet defined from an explicit set of geometry 
factors and inputs. 

• Primary inlet performance measures include: 

o Flow rates 
o Total pressure recovery 
o Inlet drag 

• Inlet aerodynamics also includes total pressure distortion, stability (buzz), etc. . . 
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Baseline Inlet Shapes 



Axisymmetric 
Pitot (K typ = 1 ) 


Two-Dimensional, Single 
Duct (K^ = 2) 



Axisymmetric, Outward- 
Turning (K^ = 3) 


Two-Dimensional, Bifurcated 
Duct (K^ = 4) 


Parts of the Inlet Flow Field 




SD Start of the subsonic diffuser 

2 Engine face 
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Freestream, Approach Flow, Engine Face 


• Freestream conditions are known from aircraft mission (Mach number, altitude). 

• The Approach Flow involves state changes (compression, expansion) due to flow over the 
forward part of the aircraft. Station L has uniform flow just upstream of the start of the inlet. 

• The Engine Face geometry and conditions (Mach number, corrected flow rate) are known as part 
of the aircraft and propulsion system design process. 


• Engine is a turbo-fan engine that typically operates at a constant corrected flow, 
approaching Mach 1 .8, the engine may reach limits that reduce the corrected flow rate. 


For M l 


Freestream 
& Approach 
Flow 
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D 2 and 
M 2 or W C2 


Engine 

Face 


External Supersonic Diffuser 


• An External Supersonic Diffuser is used for two-dimensional {K typ = 2 or 4) and axisymmetric 
outward-turning (K typ = 3) inlets. 

• The Nose is the leading edge (K typ = 2 or 4) or point ( K typ = 3) of the external supersonic diffuser. 
The nose can be sharp or round with a specified radius. 


The external supersonic diffuser consists of one or more stages through which compression 
occurs through each stage by means of an attached shock wave or isentropic Mach waves. 


Shock and Mach waves pass through focal points (e.g. shock-on-lip) 

Design strategies emphasize minimizing total pressure loss 
through the shocks and minimizing diffuser length while producing 
the desired outflow Mach number (M EX ). 

Two-dimensional diffusers include a sidewall to contain the 
compression. 


Sidewall 


Focal 

Point 


Shock 

Waves 



External Supersonic 
Diffuser 



Three-stage diffuser with shock waves 


Three-stage diffuser with smooth turning 


External Supersonic Diffuser Types 



The factor K exd controls the number of stages and design 
methods. 

K exd = 1 & 2. Explicit specification of diffuser geometry 
using A9, Ax, and Ay. 

K exd = 3. Explicit specification of diffuser geometry 
using lines, NURBS curves, polynomials, or cubic 
splines. 

K exd = 10. Single-stage ramp (2D) or cone (axi). 

K exd = 11 . Two-stage ramp (2D) or bicone (axi). 

K exd = 12. Three-stage ramps (2D only). 


Kexd ~ 

Pt 2 / Pto = 0-934 


External 

Supersonic 

Diffuser 


K exd = 1 3. Three-stages with second stage an isentropic 
compression surface. 

For 2D ramps, a carpet search method is applied to 
determine the ramp angles that result in the maximum total 
pressure recovery for the diffuser (Oswatitsch condition). 

Conical flow solutions are determined using a numerical 
solution to the Taylor-Maccoll equations. 


Bi-conic flow and isentropic diffuser contours 
determined using method-of-characteristics. 


are 


Start and end of stages are determined by path of waves 
and placement of the focal points. 



Kexd = 1 3 
Ptz / Pto = 0-964 
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Relaxed External Compression 


• Start and end of stages are determined by path of waves 
and placement of the focal points. 

• May be desirable to alter focal points and wave angles to 
provide greater margins of stability or tailor the flow. 

• Two ways of relaxing: 

1 . Distribute the focal points for the waves. 

2. Relax the angle of a downstream stage using F rejax . 
A value of F re , ax = 0 would make downstream stage 
the same angle as upstream stage. 


Individually Distributed 
Focal Points 



Waves on the cowl lip 

F relax = 0-0 

0 1 = 15° 03 = 27.3° 



Relaxed waves 

F relax = °- 5 
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Cowl Lip 


• Cowl Lip coordinate y clip is specified or determined from 
sizing operations. The code sets x dip = 0. 

• Cowl Lip planar profile can be sharp, circular, or elliptical. 



^clip — 1 

Sharp 




Kclip - 2 

Circular 








Cowl Exterior 


• Cowl Exterior provides an external surface for calculation of the wave drag. 


• Geometry model creates baseline cowl exteriors based on factor F rcex , which 
places the cowl exterior a multiple of the engine face radius. 

• Axisymmetric cowl exteriors involve extruding a profile about the axis. 



Throat 


• Throat defines interior duct from stations 1 to SD. 

• Station TH is the official “throat” station. 

• The throat cross-sections follow the shape at station 1 
(rectangular or axisymmetric). 

• Throat area distribution is set by area ratios A TH / A 1 and 
^sd / A th - 

• For K thrt = 2, the area ratio is calculated such that Mach 
number M TH is established at station TH. 


Cowl 




♦ 


(*O0cWl = (*>>’)cwlin 




External 

Si 

Di 





Subsonic Diffuser 


• The subsonic diffuser transitions from the cross-section at 
station SD to the cross-section at the engine face station 2. 

• Two-dimensional inlets require transition from a rectangular 
to a circular. Surfaces are created by the blending of 
super-ellipses. 

• Axisymmetric inlets are co-annular. Surfaces are created 
by extruding profiles about the axis-of-symmetry. 

• Geometric factors include key surface angles. 

• Length of subsonic diffuser can be directly specified or 
computed based on criteria of equivalent conical angle. 



Inlet Flow Rates 



• Quantify W h W 2 , W spillage , W bleed , W bypass , etc... 

• We know some flow rate information for the engine ( such as corrected flow, W C2 ). 

• Flow rates should satisfy continuity through the inlet streamtube. 


W, 


= 1 


W 

spillage 


+ ■ 


w 

yy i 


bleed 




+ ■ 


bypass 


W W W 

cap cap cap 


W cap is the theoretical capture flow rate. 


W r 


+ 


cap 


W 

rr jets 
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Inlet Flow Ratio = 


3 . 

w. 


(W 1 is the flow into the inlet) 


cap 


Three modes for external compression inletf low rates in supersonic flow: 

o Subcritical (inlet flow ratio < 1 , normal shock ahead of cowl lip) 

o Critical (inlet flow ratio = 1 , normal shock at the cowl lip) 

o Supercritical (inlet flow ratio = 1 , normal shock within the internal duct) 


Spillage is defined as 


Engine Flow Ratio = 


^ 7 spillage 

w. 


cap 


W r 


cap 


W, 


w: 


cap 





IP AC Manual. Fig. 3. 
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Inlet Total Pressure Recovery 


Pa 

• Inlet Total Pressure Recovery = 

PtL 

• Recovery is computed as increments through the inlet, 

Ptl _ PtEX Pr.VS Prl PtTh PrSD Pt2 
Pa. Pa Pax PrNS Pa Pan Pzsd 

• Total pressure losses modeled within SUPIN: 

o Normal, oblique, and conical shock waves 

o Viscous diffusion supersonically or subsonically 

o Sharp lip losses at high inlet flow rates at low 
speeds (high inlet flow ratio) 

• Inlet design objective is to maximize the total pressure 
recovery. 

• Characteristic cane curve relates the total pressure 
recovery to the inlet flow ratio or engine flow ratio. 



buzz limit 


subcritical 



Inlet Tools 
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Inlet Drag 


• Inlet drag consists of those drag components that can vary with inlet flow rate. 

• Inlet drag components: 

• Spillage Drag ( Additive drag + Cowl lip suction) 

• Cowl Drag ( Cowl lip drag + Wave drag) 

• Bleed and bypass drag 

• Additive drag is the sum of axial pressure forces acting on the streamtube ahead of the cowl lip. 

• Cowl lip suction accounts for positive force due to acceleration of flow around a blunt cowl lip. 

• Wave drag is the sum of axial pressure forces on the cowl exterior caused by supersonic flow. 

• Bleed and bypass drag are the result of the loss of momentum of airflow through bleed and 
bypass systems. 

• Analytical and empirical methods are used to estimate drag components. 

_ A 

^fcex 



Inlet Sizing 


• Inlet sizing establishes the inlet dimensions 
(A?ap’ I-sd> L exd , A-,, A th , A sd , angles, etc...) 

• Sizing enforces the flow rate balance through the 
inlet, and so, establishes the actual flow rates 
(W 15 W 2 ). 

• Sizing uses the total pressure recovery, which is 
dependent on the sizing of the inlet. Thus an 
iteration is involved that starts with an estimate 
of the recovery that is improved with successive 
iterations. 

• Within SUPIN, the iteration is performed by 
successive executions of SUPIN. With each 
iteration, the input data file is updated with the 
geometry properties and the total pressure 
recovery. 
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Inlet Design Study 


Apply SUPIN to design a set of baseline inlets. 

M 0 = 1 .8, h alt = 45000 ft., M 2 = 0.5, D 2 = 3.5 ft. 

Sharp nose. Three stages for the external supersonic diffuser (K exd = 13). 
M ex = 1 .3. Shock-on-lip. Sharp cowl lip. Equivalent conical angle for the 
subsonic diffuser = 3°. 


Inlet 

W 2 (slug/s) 

Pt2/Pt0 

hdip (ft) 

1- (ft) 

A C ap (ft ) 

Afcex (ft ) 

Cpwave 

Pitot 

6.4330 

0.8088 

3.196 

3.626 

8.023 

5.832 

0.3313 

2D 

7.6450 

0.9611 

2.724 

9.545 

9.534 

2.191 

0.0947 * 

Axisym 

7.6665 

0.9638 

3.489 

8.240 

9.561 

4.293 

0.2233 * 

Bifurcated 

7.6455 

0.9612 

2.724 

7.507 

9.535 

3.940 

0.0852 « 


Conclusions: 

o Pitot inlet does not perform well. 

o Flow rate (W 2 ) increases with higher total pressure recovery, but 
essentially the same between 2D and axisymmetric external compression. 

o Total pressure recovery is the same between 2D and axisymmetric 
external compression. 

o Two-dimensional bifurcated inlet is shortest inlet, 
o Cowl external wave drag is significantly less for two-dimensional inlets. 
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CFD Analyses 


• SUPIN can create planar, structured, multi-block grid for 2D or 
axisymmetric CFD. 

• Use CFD methods (Wind-US) to solve planar flow field. 

• Compute performance measures as a verification for SUPIN. 

• 3D surface grids can be start to 3D grid generation. 

Axisymmetric, Pitot Inlet (K typ = 1 ) Axisymmetric, O-T Inlet (K typ = 3 ) 


Method W 2 (slug/s) 

Pt2/Pt0 

Cpwave 

SUPIN 

6.4330 

0.8088 

0.3313 

CFD 

6.4039 

0.8094 

0.2764 

Diff 

0.45% 

-0.08% 

16.55% 


Method 

W 2 (slug/s) 

Pt2/Pt0 

Cpwave 

SUPIN 

7.6665 

0.9638 

0.2233 

CFD 

7.2946 

0.9317 

0.1706 

Diff 

4.85% 

3.34% 

23.58% 


Planar Grid 




Future Plans 


• Initial releases (March 201 1 ): 

- SUPIN as a beta test code 

- Inlet Tools eRoom 

Cruise Configuration 

- Inlet Design and Analysis Document (review draft) 


• Continue the development of SUPIN: 

- Enhance off-design performance analysis capability {K mode = 3). 

- Develop capability for design and analysis of internal and mixed- 
compression inlets. 

- Incorporate variable-geometry elements. 


Subsonic Configuration 


- Develop capability for design and analysis of three-dimensional or stream- 
traced inlets. 

- Incorporate quasi-one-dimensional CFD methods. 

• Explore other programming options and platforms (Java). 

• Explore partnerships with other researchers and organizations to 
document and develop inlet tools. 




Contact: John.Slater@nasa.aov 


(216) 433-8513 



